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Oh! Abstract: 

In exclusive J / if) production by polarized photons incident on polarized pro- 
• tons, a finite polarization asymmetry arises because of cc Fermi-motion and 

binding-energy effects. The asymmetry depends on the polarized nonfor- 
ward gluon distribution of the proton and thus gives information on gluon 
polarization in the proton. The analyzing power, however, is rather small. 
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Gluon polarization in the proton is a topic of current interest in hadronic 
physics (see e.g. |], |J). A useful probe of gluon distributions is the production 
of open or hidden charm, because charm quarks dominantly couple to gluons 
and not to light quarks in the proton. In this paper we discuss what the 
exclusive photoproduction of J/ip mesons could tell about the polarized gluon 
content of the proton. 

In general, hard exclusive processes probe nonforward matrix elements 
of the target nucleon, i.e. matrix elements of QCD operators between nu- 
cleon states of different momenta ||, These matrix elements can be ex- 
pressed in terms of nonforward parton distributions, which are functions of 
two momentum-fraction variables and generalize the usual (forward) parton 
distributions. Nonforward distributions reduce to forward distributions in the 
limit of equal proton momenta in the initial and final states. This limit can- 
not be realized in exclusive reactions because of simple kinematical reasons, 
but may be a good approximation at high energy || |||, (7|. Thus exclusive 
J/i/j production in electron-proton collider experiments, e.g. at HERA, could 
yield information on the small x behaviour of the usual gluon distribution 
of the proton, whereas fixed-target experiments will probe the nonforward 
parton distributions in a more general case. 

The dependence of unpolarized exclusive J/ip cross sections on unpolar- 
ized gluon distributions has been discussed in Refs. || §. In principle, J/xjj 
production by polarized beams on polarized targets should depend on polar- 
ized gluon distributions AG. Unfortunately, contrary to earlier belief, the 
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production amplitude does not depend on AG in the first approximation (i.e. 
assuming a nonrelativistic J/ip bound state), as we showed in a recent paper 
0- 

However, as will be shown in detail in this paper, relativistic corrections 
to J /if) production will depend on AG. Gluon polarization can therefore be 
accessed in J/if> production, although the analyzing power may be rather 
small. 

The expansion of J/if> amplitudes around the nonrelativistic limit was 
first discussed by Keung and Muzinich fTI |. Their approach is based on 
expanding the perturbative amplitudes for cc production in powers of the 
heavy quark relative momentum. Higher Fock states of the J /if} are ne- 
glected, and the calculation is therefore not gauge invariant. The issue of 



gauge invariance was discussed more recently by Khan and Hoodbhoy [12 
They pointed out that the gluonic contributions necessary to restore gauge 
inariance are proportional to v 3 , where v is the relative velocity of the charm 
quark and antiquark. Gluonic contributions are thus subleading as compared 
to the first relativistic corrections in the Keung-Muzinich approach, which 
are proportional to v 2 . 

We shall now calculate the 0(v 2 ) corrections to exclusive J /if} production 
following flTT| , |T2|| . We work in the photoproduction limit (photon virtuality 
q 2 = —Q 2 = 0) and consider only the case of collinear scattering. 

In terms of nonforward parton distributions G and AG, the helicity am- 
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plitude for 

7 (g, A) + p( Pl , S = ±1/2) - J/^K, A' = A) + p(p 2 , 5' = ±1/2) (1) 



reads 10 



1 f 1 1 

-4aa±± = - / du 



2J-i {u-i + ie){u + £-ie) 
x £ A AAAlAl £; /i 2 ) ± Ai AG(u, £; /i 2 )] , (2) 

Ai 

where £ = M 2 /{2s — M 2 ) with s = (q , ±pi) 2 . The helicity amplitude A AA / AlA2 
for 

7 (g, A) + g{h, X 1 ) -> J/^(^, A') + <?(A; 2 , A 2 ) , (3) 

where 

is the convolution of a perturbative cc production amplitude?] if AAlAa and a 
matrix element (both are Dirac matrices): 

^AA'AxAa (?) K, h) 

= Tr / 7^*^^^^' ^ ^ / rf*a:e«- (JST, A' |^(a;/2)^(-a;/2)| 0> (5) 

There are six terms in H X x 1 x 2 , corresponding to permutations of the photon 
and the two gluons on the charm quark line. A representative one, corre- 
sponding to the diagram in Fig. |T], is 

H X x 1 \ a {q,K, h'J) 

\)S F (K/2 - q + e)</(h, X 1 )S F (-K/2 -k 2 + E)f(k 2 , A 2 ) , (6) 



1 H\\ 1 \ 2 (q,K,ki;£) denotes the 73 — > ccg amplitude with heavy-quark spinors trun- 
cated, whereas in Ref. we denoted by H^,(q,K,ki) the 75 — > J/ipg amplitude with 
gluon polarization vectors truncated. 
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where Sf{p) = (J>— tn c )~ l . We have omitted normalization factors which can- 
cel in a polarization asymmetry calculation. Expanding the hard amplitude 
to second order in £, we obtain 



A AA / AlA2 (g, K, ki 
= Tr 



H XXlX2 (q,K,k i; e)\ £=0 (K, X' 



ip(x/2)ip(-x/2) 



x=0 



d 

1 d 2 

+ 2dFdi^ HxXlX ^ q,K,kl] 



x ( K, X' 



K, X' 



i/)(x/2)(-i d«)i>(-x/2) 



tfj{x/2) l-i da) [-i dpU{-x/2) 



x=0 



x=0 



(7) 



The J I ifj matrix elements can be expressed as [E^, [13 



K, X' i){x/2)i){-x/2) 

2 \ M 2 
K, X' 



K, X' 



x=0 



M 
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M 2 



tl>(x/2)[-i d a )ip(-x/2) 



x=0 



--M 3//2 
3 

4,{x/2) 



c */3 



M2 e- \g a p + ie aPlxvl » lh —j 



(9) 



6 M 2 



K a K p 



9a(3 



/* 1 + 







x=0 



(10) 



m 2 y r v" ' m / 

The amplitude given by eqs. (|7|-|T0|) is gauge invariant at 0{y 2 ). In eqs. (|8|- 
|ToD , <t> and V 2 are the J/ip meson's Coulomb-gauge wavefunction and its 
Laplacian at the origin. For a wavefunction <\> ~ e _r / a , corresponding to a 



1/r heavy-quark potential appropriate for a heavy quarkonium state, V 2 




Figure 1: One of 6 Feynman diagrams which contribute to the amplitude (|6|). 
From Ref. 

is actually infinite at the origin. The quantity V 2 in eqs. (|5|-|iT]|) has to 
be understood as representing the average of the Laplacian over a region of 
volume ~ l/m z c and thus regarded as a free parameter. 

Further corrections of the same magnitude arise because the quark mass 
m c which appears in perturbative propagators is not exactly one half of 
the physical charmonium mass which appears in the matrix elements (§-|9|). 
These binding-energy corrections are obtained by expanding the amplitude 
to first order in e# = 2m c — M. 

The parton-level amplitude simplifies to 



AxA'AxA^*?,^, fa) 
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+ 



2M M 2 J s(s - M 2 ) 



M 2 (e 1 -e*)(e,-el) 



- s( ei • e ;)(e 7 • 4) +(s- M 2 )( £l • e 7 )(e; • e 2 *)] , (11) 

where s = (q + ki) 2 and ei,e 2 ,e 7 ,e^ are the usual transverse polarization 
vectors for the gluons, the photon and the J/ip. Inserting these into the 
hadron-level amplitude (0) gives 

Aa ±± ~ f du{( ±— ^—) G ( u ^.^)[i + 0( V )} 

J-i [\u — t,+ie u + ^ — ieJ 

±A ^ il T7 - J AG( M ,e;/x 2 )) ,(12) 

2 v( M -^ + 2e ) (^ + ^-«e) 2 y J 

where 

^=M- 2 M^- (13) 
In order to estimate the magnitude of the asymmetry by using model dis- 
tributions, we shall employ the formalism of nonforward double distributions 
G(x,y) introduced by Radyushkin Q. Using 

G(u, = C dx f 1 X dy G(x, y)5{u - [x + (x + 2y - 1)£1) (14) 
Jo Jo 

and a similar relation for the polarized distribution, we can express the beam- 
target polarization asymmetry as 

da(^)-da(U) Re(P G I AG ) 

da(^)+da(U) V \Ig\ 2 1 } 



where 



Re I a = 2 f dx C * dy\n \(x + 2y) 2 - (uoxf\ dG< ^' ^ , (16) 
Jo Jo ay 



Re Jag 
lmI G 
ha- 1 AG 



rl rl—x 

— dx dy In 
Jo JO 

r2/(l+0) 

-27T / cte G(x, y)\ y 



x + 2y — ux 



d 2 AG(x,y) 



TT 





2/(l+£D) 



X + 2?/ + Co>2 

(d>-l)z/2 ' 



dx 



dAG{x,y) 



o 



<9y 



(17) 
(18) 
(19) 



y=(u-l)x/2 

(we used the notation u = !/£,). Thus the asymmetry is proportional to an 
unknown parameter 7] and depends on integrals of both G and AG. Deriva- 
tives of the distributions appear because we integrated by parts. 

The values of the parameters e#/M and V 2 0/(M 2 0) were estimated in 
where it was reported that the choice eg/M = —0.076 and V 2 0/(M 2 0) = 



—0.073 gives agreement with data on charmonium decays and inelastic J/ip 
photoproduction. In accordance with this, we shall use rj = 0.07 below. 

We now evaluate the integrals (|TB] - |T9"|) for simple model distributions and 
plot the resulting polarization asymmetry. Following the discussion of Ref. 
[ IHH , we choose 

30 



G(x,y) 
AG(x,y) 



[1-x 
30 
1 — x 



- [y(l - x - y)} xg(x) , 

- [y{\ -x- y)} 2 xAg(x) . 



(20) 
(21) 



We use the GRV-LO unpolarized gluon distribution g(x) |L5| and the Gehr- 
mann-Stirling A(LO) polarized distribution Ag(x) [p^| , evaluated at Q 2 = 
4 (GeV) 2 . Fig. g shows the asymmetry ?7Re(/^/AG)/|^G| 2 fc> r V = 0.07 as 
a function of photon energy in the proton rest frame. The asymmetry is 
of O(10~ 2 ) in the photon energy range relevant for fixed-target energies. 
At photon energies relevant for e.g. the HERA collider experiments, the 
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Figure 2: Polarization asymmetry obtained by using the model distributions of 
eqs. (|2"0"|-|2"I|) and setting 77 = 0.07, plotted as a function of photon energy in the 
proton rest frame. 



asymmetry drops by an order of magnitude. 

In summary, we have shown that a small but finite polarization asym- 
metry arises from charm quark Fermi-motion and binding-energy corrections 
to polarized exclusive J/ip photoproduction. The asymmetry depends on 
particular integrals of the polarized and unpolarized nonforward parton dis- 
tributions of the proton. Using a simple model for parton distributions and 
previous estimates of J / ip parameters, we estimate that the asymmetry may 
be of O(10~ 2 ) at fixed-target energies and O(10~ 3 ) at collider energies. 
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